A=[0 1 1 0 0;0 0 1 1 0;0 0 0 0 1;0 0 0 0 0;1 0 0 1 0];
A2=A*A;
A3=A2*A;
A4=A3*A;
A5=A4*A;
p=A+A2+A3+A4+A5;
p1=[1 1 1 1 1;1 1 1 1 1;1 1 1 1 1;0 0 0 1 0;1 1 1 1 1];
p2=p1';
p3=p1*p2;
C = intersect(p1,p2);  